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GENERAL  INTRODUCTION 

The  development  of  an  accurate  predictive  model  for  the  motion  and  deformation  of  arctic  pack 
ice  is  believed  by  most  engineers  and  scientists  concerned  with  arctic  research  and  development 
to  be  of  prime  importance  both  for  practical  operational  purposes  and  for  climatological  considera- 
tions. In  order  to  develop  such  a model,  it  was  realized  in  the  latter  part  of  the  1960’s  that  simul- 
taneous measurements  of  the  drift  of  several  different  points  on  the  ice  cover,  together  with  ap- 
propriate meteorological  and  oceanographic  measurements,  were  recpiired  (Untersteinei  and  Hnnkins 

9*).  The  basic  reason  for  this  is  that  accurate  drift  calcul  .tions  depend  upon  understanding  the 
magic tude  of  forces  caused  by  the  pack  ice  interacting  with  itself.  To  obtain  an  understanding  of 
these  forces,  it  is  necessary  to  measure  the  way  a given  area  of  ice  deforms  or  •'strains”  in  time. 
Such  defoliation  is  a result  of  differences  between  ice  velocities  at  different  lo  -ations  and  hence 
in  general  can  be  referred  to  as  differential  drift. 

One  of  the  basic  experimental  aspects  of  A1DJEX  (Arctic  Ice  Dynamics  Joint  Experiment  (see, 
or  example,  Maykut  et  al.  1972+).  has  been  to  measure  such  differential  drift  as  well  as  absolute 
drift.  In  the  springs  of  1971  and  1972.  pilot  studies  were  conducted,  one  purposes  of  which  was 
to  determine  the  optimal  way  of  sampling  the  ice  velocity  field  in  orderto  measure  the  differential 
drift  in  an  adequate  manner  for  determining  a reasonable  continuum  constitutive  law  describing 
pack  ice.  r King  these  pilot  studies  CREEL  carried  out.  among  other  studies,  a mesoscale  strain 
program.  In  this  program  the  ice  velocity  field  was  sampled  in  a region  approximately  25  km  in 
Mmeter.  The  measurements  of  the  1971  experiment  were  sparsely  spaced  in  both  time  and  space 
but  did  yield  valuable  quantitative  information  (Hibler  et  al.  1973**)  on  mesoscale  strain  not  there- 
to ore  available.  In  the  1972  experiment  a much  more  extensive  set  of  measurements  - dense  m 
oth  time  and  space  - was  obtained  which  made  possible  the  calculation  of  various  deformation 
time  series  at  3-hour  intervals  over  a period  of  more  than  30  days.  The  insults  allowed  a detailed 
study  to  be  made  of  the  nature  of  the  ice  velocity  field.  Also,  using  available  meteorological  data 
the  results  were  compared  with  theoretical  ice  drift  models  with  good  success. 

This  report,  divided  into  two  parts,  describes  the  analysis  of  the  1972  results.  Fart  1 examines 
the  results  with  a view  toward  understanding  of  the  nature  of  the  ice  velocity  field  and  determining 
t .3  minimum  scale  on  which  the  pack  ice  may  be  viewed  as  a continuum.  Part  II  represents  the 
formulation  of  an  analytic  solution  of  a linear  drift  model  which  includes  ocean,  air  and  Coriolis 
forces  and  models  the  ice  as  a viscoelastic  medium.  The  predicted  results  of  this  calculation 
ased  on  meteorological  variables)  are  compared  with  observed  differential  drift  measurements  with 

good  agreement  and  a number  of  conclusions  regarding  ice  drift  calculations  using  a linear  model 
are  drawn. 
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**  Illblfir-  W n ,11-  WF.  Weeks.  S.  Ackley,  A.  Kovacs  and  W.J.  Campbell  (1973)  Mesoscale  strain  measure 
ments  on  the  Beaufort  Sea  pack  ice  (AIDJF.X  1971).  Journal  of  Glaciology,  vol.  1L_\  no.  J p.  187.-06. 
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by 

W.D.  Ilibler  III,  W.F.  Weeks,  A.  Kovaes  and  S F.  Ackley 

Introduction 

One  of  the  primary  goals  of  the  Arctic  Ice  Dynamics  Joint  Experiment  (AIDJKX)  is  an  improved 
understanding  of  the  drift  of  pack  ice.  To  tins  end  accurate  field  observations  of  the  deformation 
of  different  types  of  paek  ice  performed  on  a variety  of  time  and  space  scales  are  urgently  needed. 
To  pat t tally  satisfy  this  need  a series  of  detailed  mesoseale  strain  measurements  were  made  at 
approximately  3-hour  intervals  over  a 30-day  period  in  the  spring  of  1972  in  the  Beaufort  Sea.  These 
measurements  are  particularly  useful,  since  in  outlier  studies  of  sea  ice  deformation  - as  reviewed, 
for  example,  by  Hiblor  et  al.  (1973a)  - there  were  usually  large  and  inudom  time  intervals  between 
observations,  making  the  computation  of  accurate  time  series  impossible.  Also,  and  perhaps  even 
more  important,  thesr  earlier  studies  included  no  detailed  investigation  of  the  nonlinear  variations 
in  the  ice  velocity  field  that  result  from  inhomogeneous  spatial  variations  or  fluctuations  in  the 
deformation  of  the  ice. 

Therefore,  the  analysis  described  here  was  undertaken  with  three  primary  goals  in  mind  first, 
to  provide  a detailed  time  scries  of  the  least  squares  strain  rate  tensor  and  vorticity  (complete 
with  "error  bars"  due  to  the  nonlinearity  of  the  ice  velocity  field)  over  the  25-dav  period,  Julian 
day  88-113  (28  Mar-22  Apr).  1972,  second,  to  study  the  magnitude  and  nature  of  the  nonlinear  ve- 
locity fluctuations;  and  third,  to  compare  deformation  measures  from  different  scales  to  determine 
coherent  modes  of  deformation  in  different  sized  arrays  as  well  as  sealing  effects.  These  results 
can  then  be  compared  with  predictions  from  theoretical  drift  calculations  and  with  data  collected 
on  the  remote  sensing  overflights.  Besides  providing  insight  into  the  nature  of  pack  ice  dynamics, 
such  comparisons  and  information  on  inhomogeneities  in  the  ice  velocity  field  are  helpful  in  de- 
signing future  strain  arrays. 


Approach 

As  a framework  in  which  to  view  our  analysis,  it  is  useful  to  think  of  the  paek  as  a large  num- 
ber of  irregular  ice  floes  packed  closely  together,  with  the  compactness  varying  with  season.  In 
the  summer,  when  the  compactness  is  low,  the  individual  ice  floes  can  readily  be  identified  and  the 
pack  looks  like  a two-dimensional  granular  medium.  In  the  winter,  when  the  compactness  is  near  one, 
the  individual  floes  can  no  longer  lie  clearly  identified  and  the  ice  is  crisscrossed  by  a number  of 
irregular  leads.  A typical  example  of  the  pack  ice  structure  in  winter  is  given  in  Figure  1,  which 
shows  a schematic  diagram  of  active  leads  and  ridging  zones  in  the  mesoseale  stnin  array  for  one 
instant  of  time. 

Using  such  a conceptual  model  of  pack  ice,  the  velocity  of  any  point  in  the  paek  can  be  viewed 
as  consisting  of  a continuum  velocity  component  (varying  over  lengths  commensurate  with  the  scale 
of  meteorological  variations)  plus  a fluctuation  component  due  to  the  discrete  small-scale  nature  of 
the  paek  ice.  Such  a partition  of  the  velocity  field  is  similar  to  that  which  can  be  made  for  a fluid 
where  the  motion  of  each  molecule  has  a continuum  component  plus  a fluctuation  component.  In  the 
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Figure  1.  Schematic  diagram  of  the  mesoscale  strain  array  together  with  active  leads  and  ridging 
zones.  Leads  and  ridges  were  obtained  from  a 1500-m  "rial  photo  mosaic  taken  on  6 April  1972. 

case  of  the  pack  ice,  however,  the  fluctuations  are  probably  m due  so  much  to  random  motion  of 
floes  (although  this  is  undoubtedly  a factor)  as  to  the  fact  that  the  floe  sizes  are  large  relative  to 
our  measurement  scale  and  thus  spatial  velocity  profiles  are  of  a stepwise  nature.  This  is  illus- 
trated by  Figure  1,  which  shows  that  the  distance  between  leads,  whore  relative  motion  occurs,  is 
generally  4 to  6 km.  Also,  for  pack  ice,  the  fluctuations  are  expected  to  be  highly  variable  in  time, 
since  they  are  primarily  driven  by  the  transfer  of  energy  into  the  pack  by  meteorological  forces. 

This  is  opposed  to  a fluid  in  a laboratory  in  which  fluctuations  are  well  described  by  the  tempera- 
ture. 

As  we  examine  spatial  variations  in  the  velocity  field  over  larger  and  larger  areas,  the  con- 
tributions of  fluctuations  to  average  velocity  differences  (i.e.  strain  rates)  would  be  expected  to  be- 
come less  pronounced.  Stated  differently,  the  contribution  of  fluctuations  to  the  sea  ice  strain  rate 
would  be  expected  to  become  small  when  the  area  covered  by  the  strain  array  becomes  large  relative 
to  floe  size  and/or  distances  between  leads. 

In  order  to  sort  out  the  fluctuations  from  the  continuum  motion,  we  analyzed  the  position  data 
of  the  mesoscale  targets  by  fitting  a least  squares  planar  surface  to  the  spatial  velocity  field 
sampled  by  the  array,  with  the  slope  angles  of  the  plane  representing  the  strain  rate  and  vorticity. 
Such  a procedure  is  commensurate  with  the  discussion  of  sea  ice  strain  by  Nye  (1973),  who  notes 
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th.it  tor  a pice  1 sc  measure  of  the  strain  rate  tin*  velocity  field  needs  first  to  he  smoothed  hi  foie 
d nvatives  are  taken.  Since  the  area  covered  by  the  niesoscale  array  is  relatively  small  ( 20  km 
in  diametei ) com])  ired  with  inoteoiologi'al  systems,  we  would  expect  the  continuum  \ i locity  tohe 
iclatively  linear  over  this  legion  and  hence  the  planar  approximation  should  he  (rood.  Residual 
(h  \ i at  ions  tiom  this  planar  surface  arc  identified  as  fluctuations.  These  residual  deviations  also 
cause  some  uncertainty  in  the  strain  rates,  an  uncertainty  which  can  lx;  estimated  by  dimensional 
analysis  for  different  sized  ariays  on  the  assumption  that  the  fluctuation  amplitude  is  reasonably 
similai  over  the  legion.  As  part  of  this  proeeduie  to  delineate  the  c.fference  between  continuum 
motion  and  fluctuations,  we  also  estimated  a continuum  length,  defined  as  the  length  ovei  which, 
on  the  average,  the  continuum  velocity  differences  equal  the  fluctuation  amplitude.  This  length, 
which  is  about  10  km,  gives  a lough  measuie  of  the  scale  above  which  sea  ice  may  he  viewed  as 
a continuum  and  below  which  the  diseiete  nature  of  the  pack  begins  to  dominate  the  motion. 

Finally,  in  addition  to  such  a h ast  squares  analysis  for  determining  scaling  effects,  we  also 
1)  computed  strains  obtained  from  triangles  of  diffeient  sizes  ( 5 km  to  20  km)  and  2)  computed 
the  niesoscale  strain  lesults  over  the  20-km  region  with  the  macroscale  strain  results  obtained  fiom 
a 100-km  triangle,  one  corner  of  which  was  the  cento  of  the  niesoscale  array.  These  comparisons 
genet  ally  indicated  that  all  arrays  were  measuring  similar  continuum  motions  of  the  pack,  with  the 
fluctuations  yielding  a huge  contribution  to.  but  not  masking,  the  continuum  motion  on  a scale  of 
about  20  km. 


Site  Location  and  Data  Collection  Procedures 

The  rut  asurements  used  for  this  study  were  made  m the  vicinity  of  the  mam  1072  A1DJKX  camp, 
located  at  roughly  75  00'N.  MS  30 W.  The  camp,  as  well  as  the  diffeient  research  programs  con- 
ducted fiom  it.  are  desenbed  in  A1D.JKX  Bulletin  no.  11  (Thorndike  et  al.  1972).  The  strain  array 
was  established  by  erecting  a series  of  targets  which  consisted  of  eoinei  cubes  mounted  on  the  tops 
of  aluminum  poles.  The  distanct  s to,  and  angles  between,  the  taigets  were  measured  using  a con- 
tinuous wave  laser  langefindei . The  height  of  the  targets  vaiied  from  3 to  10  m above  the  ice  sin- 
face  depending  on  ihe  distanct  and  the  obstructions  between  the  targets  and  the  main  camp.  A 
(liagiani  of  the  sttaiti  allay,  together  with  an  oveilay  of  active  leads  and  ridging  zones,  is  shown  in 
Figure  1.  The  angles  between  the  targets  weie  measured  with  ail  average  accuracy  of  bettei  than 
’ 1 min  and  were  referenced  to  a fixed  stake  on  the  multiyear  floe  on  which  the  main  camp  was 
sited.  The  line  between  the  laser  and  the  stake  was  then  tied  into  the  tine  north  dotei munitions 
(sun  shots)  made  by  Thorndike  and  (lill  (Thorndike  et  al.  1972)  Distances  were  measured  to  the 
neatest  0.03  m (0.1  ft)  because  the  huge  strains  that  were  experienced  obviated  the  need  for  any 
gteater  piecisioti. 

This  stiain  measurement  system  was  found  to  be  vastly  superior  to  manned  tclliuonieter  sites 
(lliblcr  et  al.  1973a).  With  this  system,  a large  inimlier  of  strain  lines  could  be  determined  easily 
without  manning  the  i emote  stations.  It  was  also  iclatively  quick  and  easy  to  install,  and  placed  a 
minimal  reliance  upon  "black  boxes."  llowcvet.  visibility  problems  (wind-blown  snow,  sea  smoke 
fiom  leads)  made  acquisition  of  continuous,  equally  spaced  time  si  ties  difficult  (laser  measurements 
weie  impossible  approximately  10T  of  the  time,  and  once  high  winds  caused  a gap  of  almost  two 
days  in  the  strain  line  time  series).  In  addition,  the  system  required  manpowei  21  hours  a day 
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Figure  2.  Typical  results  of  the  interpolation  and  smoothing  process  used  to  generate  ei/uispaced 
values  lor  the  strain  analysis.  Curve  a was  obtained  with  a smoothing  filter  transition  band  from 
S.O  to  6.15  hours,  and  curve  b with  a transition  hand  from  20  to  11  4 hours. 

Data  Analysis 

Extrapolation  and  smoothing  of  strain  data 

Data  taken  in  the  field  consisted  of  distances  and  angles  of  targets  relative  to  a fixed  reference 
stake.  As  a first  step  in  the  data  (eduction,  these  angles  were  converted  to  angles  relative  to  true 
north  as  measured  on  Julian  day  81  (21  March).  Rotations  of  the  array  were  not  taken  into  account 
for  strain  calculations,  so  that  the  coordinate  system  used  for  this  study  is  slightly  different  from 
the  true  north  coordinate  system  (the  maximum  difference  is,  however,  less  titan  5 ).  Rotations  of 
the  array  were  included  in  the  vorticity  calculations  The  time  scale  was  converted  to  GMT  by 
using  four  GMT  calibration  times  obtained  in  the  field  to  find  a least  squares  rate  for  the  clock  used 
in  the  measurements.  The  data  point  times  were  recomputed  with  tins  rate  and  then  all  data  (both 
angles  and  distances)  were  linearly  interpolated  and  resampled  evety  hour.  Using  tins  new  data 
set  at  1-hour  intervals,  the  time  series  was  smoothed  with  a lo  pass  filter  having  a transition  band 
width  with  periods  from  8.0  to  6.15  hours.  The  smoothed  time  series  was  then  resampled  every 
third  hour.  If  there  was  reason  to  expect  that  the  data  contained  a large  number  of  time  intervals 
greater  than  3 hours,  a low  pass  filter  having  a transition  band  from  20  to  11.4  hours  was  used  be- 
fore resampling.  Both  filters  had  less  than  0.6%  side  lobe  errors  and  consisted  of  81  symmetric 
weights  designed  according  to  the  procedure  discussed  by  Ilibler  (1972). 

This  process  of  interpolation,  followed  by  smoothing,  maybe  viewed  as  a consistent  way  of 
constructing  a smooth  curve  (with  no  high  frequency  components  past  a reasonable  cutoff  dictated 
by  the  average  sampling  rate)  through  the  randomly  spaced  data  points.  Alternatively  the  curve 
may  be  considered  an  accurate  representation  of  the  low  frequency  portion  of  the  linearlv  interpolated 
curve. 

Examples  of  the  smooth  curves  generated  by  this  process  are  shown  in  Figure  2.  Curve  a.  which 
results  from  the  filter  with  the  higher  frequency  cutoff,  follows  the  data  quite  closely.  This  indu  ales 
that  there  is  little  variance  associated  with  periods  shorter  than  6 hours. 
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Experimental  error  estimation 

The  primary  source  of  error  in  the  measurements  was  the  uncertainty  of  target  angles.  Since 
angular  measurements  were  generally  accurate  to  t 0,0003  radians  (4  1 min)  and  distance  errors 
were  small,  we  estimated  the  x and  y position  errors  of  each  measurement  of  a target  at  distance 
r and  angle  0 to  be 


\x  0,0003  r sin  6 
\ y 0,0003  r cos  0 . 


(1) 


Since  measurements  3 hours  apart  were  subtracted  to  obtain  velocities  v,  we  estimated  x and 
y velocity  errors  to  be 


Av  0.0006  r sin  0/At 

(2) 

A vy  0 0006  r cos  0/ A t 

with  \t  3 hours.  This  is  a slight  overestimation  of  errors  for  the  difference  between  two  numbers 
with  uncorrelated  errors.  However,  since  there  may  have  been  some  errors  - due  to  interpolation 
between  nonequally  spaced  points  - that  were  not  removed  by  smoothing,  we  have  made  a con- 
servative error  estimation.  The  experimental  errors  given  by  eq  2 were  then  used  as  input  to  obtain 
the  experimental  error  on  the  strain  rate  tensor,  as  discussed  in  the  next  section. 


Least  squares  comuutational  technique 

To  understand  conceptually  the  least  squares  strain  and  vorticity  calculations,  it  is  useful  to 
visualize  a contour  plot  of  the  x (or  y)  velocity  component  of  the  ice.  In  essence,  the  computer 
program  used  to  calculate  the  deformation  rate  fits  a planar  surface  through  the  eor.tour  plot  with 
the  slope  angles  of  the  planes  yielding  the  strain  rates  and  vorticities.  Since  the  aetual  velocity 
components  deviate  from  a perfect  plane,  there  is  some  uncertainty  in  the  slope  angles  of  the  plane. 
We  refer  to  the  average  deviation  of  the  velocity  components  from  the  plane  as  the  residual  fluctua- 
tion error  and  the  uncertainty  in  the  slope  angles  of  the  plane  as  the  inhomogeneity  variation.  In 
addition,  once  the  least  squares  equation  for  the  plane  as  a function  of,  say  JV,  velocity  measure- 
ments is  determined,  the  estimated  experimental  errors  may  be  used  as  input  to  obtain  estimates  of 
slope  uncertainties  due  only  to  experimental  errors.  The  remainder  of  the  slope  uncertainty  may 
then  be  identified  with  nonlinearities  in  the  iee  velocity  field. 

To  formulate  this  conceptual  model  mathematically,  we  proceed  as  follows.  Using  tensor  nota- 
tion for  the  strain  rate,  the  strain  rate  tensor  and  vorticity  w are  defined  by 

; 

1 f'lv2  ,lvl\ 

’ 

where  v is  the  ith  velocity  component  of  the  ice  pack  (considered  as  a continuum)  and  i,  j 1,2, 
since  we  are  only  concerned  with  the  horizontal  motion  of  the  sea  ice.  Considering  N targets  whose 
positions  are  being  measured,  we  denote  by  vj'*  the  measured  jth  velocity  component  of  the  ith  tar- 
get, and  by  0,  and  r(  the  polar  coordinates  of  the  ith  target  relative  to  an  arbitrary  origin. 
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As  a model  io  explain  the  velocities  of  the  N targets,  we  consider  the  ice  velocity  field  (on 
the  scale  of  about  20  km)  to  consist  of  a continuum  velocity  component,  which  varies  in  a uniformly 
linear  manner  (specified  by  < and  w),  and  a random  fluctuation  ^ nmponent.  Mathematically  our 
model  is  expressed  by  the  equation 


Vj 


4 Z 


where 


(3) 


v\n,  i 

1 to  N 

II 

1 

1 

vO 

N\  < 

N t 1 to  2N 

Uj  cos  0{ 

r(  sin  0 ■ 0 

0 

1 

0).  i 

1 to  N 

X.J 

[ o 

0 fj  cos  0 

rt  sin  0X 

0 

11.  i 

N ♦ 1 to  2N 

( < 
J 

' 

?ll 

;lL>  - w 

f,2  ' w 

'22 


and  are  constants  representing  the  continuum  vc’oeity  components  a!  the  origin.  In  eq  3, 

Zj  is  the  fluctuation  component  plus  any  measurement  error  and  E(ut)  (■  is  the  “expected" 
value  of  u(,  since  E(z ,)  0.  The  least  squares  estimates  of  denoted  by  <\  are  obtained  by 

* p * * 1 

minimizing  1 z‘. 

i 2N 

To  do  this,  we  differentiate  1 (u,  - f.x,,)2  with  respect  to  t.,  which  yields  the  matrix  equa- 

j 1 1 J * 

tion  for  the  least  squares  estimates  of  <)  (Jenkins  and  Watts  1969,  p.  131): 


('=M~lX'U 


(4) 


where  M X'  X and  prime  denotes  transposes. 

When  using  these  least  squares  equations,  note  that  adding  a constant  rotation  to  all  angles 
only  changes  the  vorticity  w.  This  can  be  demonstrated  by  noting  that  in  eq  3 changing  w to  w *■ 
0 changes  ut  to  u/  where 

uj  - (sin  0rf  i = 1 to  N 

ui'  B 

u(  + (cos  0j)  0r}  i N r 1 to  2N  . 
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This,  however,  is  equivalent  to  adding  a constant  rotation  to  all  angles.  In  a similar  manner,  add- 
ing a constant  velocity  to  all  points  changes  only  A , and  A„. 

Since  we  have  only  a finite  number  of  random  velocity  measirements  (each  with  some  random 
fluctuation  “error”)  to  calculate  <(,  there  is  some  uncertainty  or  error  in  r . To  estimate  this  vari- 
ation or  uncertainty,  it  is  necessary  to  calculate  the  covariance  matrix  of  <*,  which  is  easily  shown 
to  he  given  by  (Jenkins  and  Watts  1969,  p.  134) 

C ')'  (5) 

where  cov(<^,  f^)  and  co v(zjt  z-).  In  our  case  z.  consists  of  two  parts 


where  z*1  is  the  measurement  error  (which  we  estimate  to  be  given  by  eq  3)  and  z[  is  the  fluctuation 
component  of  the  ice  motion,  since  the  actual  ice  velocity  field  is  nonlinear.  If  we  assume  that 
zf  and  z[  are  uncorrelated,  then  the  variation  in  r due  only  to  measurement  errors,  call  it  C^.  may 
be  estimated  using  eq  5 with  co v(zjf. /“).  Using  eq  2 we  estimate  Fjj  to  be 


VM 

ij 


(sin  Of  r2(\0  \tf. 
* (cosf^)2  r2  (A  0 \t)2, 


i 1 to  N 

i /V  * 1 to  2/V. 


(6) 


To  find  the  variation  of  some  linear  combination  of  the  r]  due  to  the  total  “error’’  z , we  make 
the  usual  assumption  that  z ■ are  uncorrelated  with  the  same  mean  and  variance  so  that  t ov(z,,  z ) 

c 0 1 o o 1 J 

where  the  point  estimator  of  a (denoted  by  s ) is  given  by  the  residual  error 


2N  - 6 


(7) 


where  u;  =■  . 

In  this  case  eq  5 reduces  to 

Cij  =(|r1)ij*2.  (8) 

^ ZN  A zN  /\ 

Given  a linear  combination  of  f(,  h n^,  then  the  deviation  of  h from  E(b)  is  such  that 

b - E{b)/sy7^W\~  (9) 

has  the  t distribution  with  2N- 6 degrees  of  freedom,  assuming  that  the  errors  are  normally  distributed 
(Bennett  and  Franklin  1954,  p,  250).  Consequently,  confidence  limits  for  the  estimated  strain  may 
be  obtained  using  a t distribution  table. 

In  addition  to  the  estimated  strain,  we  are  also  concerned  withthe  velocity  fluctuation  compo- 
nent zj,  which  is  strictly  speaking  not  an  error,  but  represents  the  variation  from  linearity  of  the 
velocity  field  over  the  region  sampled.  In  the  cases  we  have  studied,  the  estimated  residual  error 

z,  obtained  from  eq  7 was  generally  found  to  be  larger  than  the  average  estimated  value  of  experi- 
1 M 

mental  error  z;  from  eq  2.  Consequently,  as  a matter  of  convenience,  we  w>ll  often  refer  to  the 
residual  error  obtained  from  eq  7 as  the  residual  fluctuation  error. 
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Also,  as  a matter  of  notation,  we  will  refer  to  the  uncertainty  in  (\ttovT^T)  as  the  mhonio- 
geneity  error.  It  should  he  lememliered  that  the  liihomogeneity  error  is  an  estimate  of  the  un- 
certainty in  the  least  squares  estimated  strain  and  consequently  depends  somewhat  on  the  number 
of  samples  used.  The  residual  fluctuation  error,  oil  the  other  hand,  should  be  relatively  indepen- 
dent of  the  number  of  sample  points  used. 


Strain  Results 


For  comparison,  strain  rate  and  vorticity  time  series  have  been  calculated  using  several  dif- 
ferent combinations  of  targets.  Table  1 describes  all  combinations  of  targets  used  in  the  calcula- 
tions. The  origin  was  also  considered  a position  measurement  in  some  arrays,  with  distance  and 
angle  being  zero  for  all  time. 

Over  the  time  interval  studied  in  this  paper,  there  was  only  one  major  gap  m the  time  series, 
and  this  portion  is  blanked  out  in  the  plots,  However,  lor  spectral  studies,  root  mean  square  (mis) 
error  estimation,  and  correlation  studies,  the  whole  curve  was  used  (data  points  every  3 hours', 
including  linearly  extrapolated  data  through  t lie  gap,  with  the  linear  extrapolation  being  done  on 
each  target  position  as  discussed  previously. 


Table  I.  Strain  line  combinations  used  in  this  paper. 


Array 

"name" 

Targets 

Approximate  diameter 
ol  region  sampled 
(km) 

Combined  niray 

i,  2,  a.  8,  n,  4.  7,  li, 
13,  5,  10,  origin 

20 

Oilier  array 

1,  2,  3,  8,  0,  origin 

20 

16 -km  triangle 

2.  5.  ') 

lb 

8-km  triangle 

11,  13.  4 

8 

r>-km  inangte 

11,  7,  origin 

f> 

Strain  tensor  time  series 

Since  the  outer  array  targets  were  measured  more  often  than  other  targets,  the  outer  array 
provides  a more  detailed  time  series  for  analysis,  Results  of  least  squares  calculations  using  the 
outer  array  are  shown  in  Figures  3 and  4,  which  present  the  two  invariants  (Nye  1957)  (denoted  bv 
f j and  <„)  of  the  strain  rate  tensor  both  separately  and  in  the  form  of  the  divergence  rate  fn  (sum 
of  the  principal  axis  components)  and  the  maximum  shear  rate  (<  t <o)/2  (difference  between  the 
principal  axis  components  divided  by  2).  For  the  strain  rates,  the  angles  of  the  strain  lines  at  the 
beginning  of  each  3-hour  interval  were  used  m the  least  squares  calculations.  The  continuous 
errors  in  Figures  3 and  4,  the  \ values  shown  below  each  curve,  represent  the  inhomogeneity  error 
(calculated  using  eq  7 and  8),  which  is  due  primarily  to  velocity  fluctuations.  The  small  error  bar 
on  the  divergence  rate  curve  represents  the  maximum  experimental  error,  which  was  calculated  by 
first  calculating  the  experimental  error  for  each  point  in  the  time  series  according  to  eq  5 and  6 
and  then  finding  the  maximum  error  over  the  time  series.  For  a more  compact  summary  of  the  rela- 
tive magnitude  of  the  strain  rates  and  strain  rate  variation  errors,  we  list  the  root-mean-square 
(rms)  values  for  the  various  time  series  in  Table  II. 
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Figure  3.  Least  squares  divergence  rate  and  accom- 
panying inhomogeneity  error,  and  maximum  shear  rate 
and  inhomogeneity  error.  The  small  error  bar  in  a 
represents  the  maximum  uncertainty  due  to  measure- 
ment error. 
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Figure  4.  Principal  axis  components  of  the  least 
squares  strain  rate  tensor  and  inhomogeneity  errors. 
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Table  II.  Root- mean-square  strain  rates,  strain  rate  variation 
errors  and  experimental  errors  (or  combined  array. 

(Units  10  4 hr-’). 


N,'l  •'«„  A«V  (expert  (.‘j-fJ/2.  \U , Ij  2| 

10. 'i  1.92  3.66  1.04  3.28  2.78  0.85  3.48  1 30 


From  Figures  :!  and  1 and  Table  II.  we  see  that  m general  stiain  rates  are  of  greater  magni- 
tude than  their  respective  strain  rate  inhoinogeneity  errors.  The  figures  also  indicate  that  the 
lnhoniogeneity  errors  generally  increase  with  increasing  strain  rate.  It  is  important  to  remember 
that  the  inhomogeneity  errors  shown  are  for  6 targets  and  would  tie  smaller  for  a larger  numtier  of 
targets.  This  is  analogous  to  the  error  on  the  slope  of  a simple  least  squares  line  which  becomes 
less  as  more  prints  are  added  even  though  the  standard  error  of  the  estimate  may  remain  the  same. 
In  particular,  for  the  same  residual  error  for  each  velocity,  an  equilateral  strain  triangle  would  have 
a divergence  rate  error  1.7  times  as  large  as  that  shown  in  Figures  3 and  1. 

One  striking  aspect  of  the  deformation,  best  illustrated  in  Figure  4.  is  that  m the  principal 
axis  coordinate  system  most  of  the  expansion  or  contraction  is  taking  place  along  one  axis.  More- 
over. there  is  usually  contraction  along  one  axis  and  extension  along  another.  Another  salient 
characteristic  of  the  ice  deformation  is  that  the  ice  motion  appears  to  consist  of  deformat lonal 
events  which  occur  every  several  days  and  usually  consist  of  dilation  followed  by  convergence. 

The  strain  rate  time  series  generally  shows  rather  rapid  variations  in  strain  which  aie  probably 
due  to  the  random  bumping  and  yielding  of  floes  as  well  as  the  random  opening  of  leads.  Under 
our  idealized  model,  consisting  of  ice  fluctuations  superimposed  upon  a continuum,  these  high 
frequency  motions  should  primarily  represent  fluctuations.  With  respect  to  maximum  strain  rates, 
the  maximum  observed  divergence  rate  is  0.16"'.  * 0.09",'.  per  hour  with  a maximum  convergence  rate 
of  0.15%  ♦ 0.09%  per  hour.  The  largest  maximum  shear  r ite  is  0.16%  ♦ 0.05%  per  hour. 

Nonlinear  velocity  fluctuations 

The  least  squares  calculation,  besides  yielding  the  average  strain  rates,  also  gives  a mea- 
sure of  the  nonlinearity  in  the  velocity  field  through  the  residual  fluctuation  error.  This  residual 
erroi  can  be  viewed  as  a fluctuation  in  the  velocity  field  from  the  ideal  continuum  value.  The 
magnitude  of  these  fluctuations  is  important  in  determining  the  size  of  a measurement  array  neces- 
sary to  measure  accurately  the  average  strain  rate.  In  terms  of  our  continuum  model,  the  fluctua- 
tion contribution,  together  with  the  average  strain  rate,  yields  a characteristic  length  above  which 
the  pack  ice  may  he  considered  a continuum  and  under  which  the  ice  motions  of  individual  floes 
and  leads  become  dominant.  Such  a characteristic  length  is  estimated  by  determining  the  length 
over  which  the  fluctuation  on  the  average  is  almost  the  same  as  the  continuum  velocity  variations. 

For  a best  estimation  of  the  fluctuation  error,  we  utilized  the  combined  array  consisting  of  11 
targets  plus  the  origin.  Since  some  targets  in  the  combined  array  were  not  measured  as  frequently 
as  those  in  the  outer  array  the  linearly  extrapolated  data  were  smoothed  with  a filter  having  a 
transition  band  from  20.0  to  11.4  hours,  as  discussed  in  the  previous  section  on  data  analysis 
(p.  5). 

The  resulting  residual  errors  from  the  combined  array,  and  strain  tensor  components  in  a 
(north,  west)  coordinate  system  are  shown  in  Figure  5.  In  order  to  put  the  residual  error  and  de- 
formation rates  in  perspective  with  the  overall  motion  of  the  pack  ice.  plots  of  the  velocity 
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Figure  5.  Comparison  of  mcsoscale  divergence  rate, 
north-south  strain  rale,  east-west  strain  rate,  shear 
rate,  least  squares  residual  fluctuation  error,  east- 
west  velocity  of  the  camp,  north-south  velocity  of 
the  camp,  and  speed  of  the  camp.  All  curves  were 
smoothed  with  a low  pass  filter  having  a transition 
band  from  20.0  to  11.4  hours. 
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components  of  the  central  point  of  the  array  are  also  shown  in  Figure  5.  These  velocity  plots 
were  obtained  using  drift  data  obtained  by  Thorndike  (1973)  from  satellite  navigation  fixes.  The 
velocity  data  were  smoothed  with  the  same  filter  used  in  the  strain  calculation.  In  the  processing 
of  the  data  carried  out  by  Thorndike,  some  smoothing  was  also  done,  but  this  smoothing  affected 
only  the  higher  frequencies  past  the  transition  band  of  our  filter.  For  a more  concise  comparison 
of  the  data,  rms  values  of  the  various  curves  are  given  in  Table  111. 

As  can  be  seen  from  Table  111,  the  characteristic  length  (residual  error/strain  rate)  over  which 
the  velocity  change  due  to  strain  rate  is  of  the  same  magnitude  as  the  residual  fluctuation  error  is 
of  the  order  of  10  km.  A second  characteristic  length  of  some  interest  is  the  length  over  which  the 
velocity  difference  according  to  the  strain  rate  equals  the  average  drift  velocity.  This  length  is 
of  the  order  of  1000  km,  roughly  similar  to  the  size  of  the  Pacific  Gyre.  Consequently,  the  residual 
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Table  111.  Root-aie aa-equare  strata  rates,  residual  errors, 
and  central  point  velocities  (from  Pig.  5). 


D iv  rare 

(hr'1)  • j | (hr'  *)  i^,  (hr  ') 


• ls(hr  li 


Residual 
error 
(m  hr 


East-west 

vet 

( m lir~  *) 


North-south 

vel  Speed 

(ro  hr"1/  (m  hr  *) 


i.-'kiiH)4)  2.60<io4)  2.37<to~4) 


2.12  < in-4) 


291;.  3 


107.9 


312.9 


''trot  is  rattier  insignificant  in  terms  of  the  absolute  velocity  for  a point,  but  becomes  more  critical 
in  terms  of  the  velocity  gradient  or  deformation  rate.  In  terms  of  a continuum  model,  the  results 
suggest  that  on  a s<  ale  larger  than  about  10  km  the  pack  ice  begins  to  behave  like  a continuum, 
whereas  on  a smaller  scale  the  individual  particle  behavior  begins  to  dominate  the  observed  motion. 

The  magnitude  of  the  residual  error  also  allows  an  estimate  to  be  made  of  the  effect  of  fluctua- 
tions on  strain  rates  from  arrays  of  different  sizes.  For  example,  we  can  view  the  strain  measured 
by  a simple  area,  say  a triangle,  to  consist  of  two  components:  1)  the  continuum  strain  rate  and 
'•!)  the  fluctuation  component.  Assuming  that  the  residual  error  is  approximately  constant,  inde- 
pendent of  the  size  of  the  triangle,  then  the  fluctuation  contribution  to  the  strain  rate  is  of  greater 
relative  magnitude  for  small  triangles.  In  addition,  since  the  fluctuation  component  would  be  ex- 
pected to  consist  of  rather  rapid  bumping  motions,  strain  rates  would  behave  in  a mere  erratic 
manner  as  the  measuring  triangle  becomes  smaller.  This  effect  is  apparent  in  plots  of  individual 
triangles.  Figure  6.  for  example,  illustrates  the  net  divergences  (essentially  the  areas)  of  throe 
nested  triangles  as  a function  of  time.  The  more  rapid,  large-magnitude  motion  of  the  small  tri- 
angles is  apparent.  The  curves  do,  however,  illustrate  a general  correlation  of  strain  events  con- 
sisting of  dilation  followed  by  convergence.  Generally  tiiese  results  suggest  that,  although  the 
smaller  triangies  are  averaging  over  very  few  leads  and 'or  floes,  over  a period  of  several  days  the 
ice  on  a small  scale  would  be  expected  to  diverge  if  the  pack  is  generally  diverging  and  converge 
if  the  pack  is  converging. 


Comparison  of  Mesoscale  Deformation  with  Macroscale  Deformation 

In  addition  to  a satellite  navigation  system  at  the  main  1972  AID.JF.X  camp  (the  center  point 
of  the  inesoscalc  array),  the  1972  pilot  study  included  satellite  position  measurements  of  two  other 
camps  approximately  100  km  west  and  northwest,  respectively,  of  tin*  main  camp.  Strain  data  from 
this  larger  triangle,  referred  to  as  macroscale  deformation,  provide  a valuable  measure  of  the  de- 
fixmation  averaged  over  a larger  region  than  that  covered  by  the  mesoscale  array. 

In  a comparison  of  the  macroscale  and  mesoscale  deformation  rates,  arguments  for  both  simi- 
larities and  differences  can  be  made.  First,  since  weather  systems  typically  vary  over  several 
hundred  kilometers,  we  would  intuitively  expect  some  similarity  between  the  macro-  and  mesoscale 
deformation  rates,  at  least  to  the  extent  that  both  systems  an-  measuring  the  continuum  motion  of 
the  ice  pack.  However,  there  are  several  reasons  why  the  correlation  should  not  be  extremely  good 
Foremost  is  the  fact  that  the  mesoscale  array  is  only  slightly  larger  than  the  estimated  ontinuum 
length  of  10  km  and  tints  fluctuations  of  tin*  velocity  field  yield  a strong  component  in  the  least 
squares  mesoscale  deformation  time  series,  a much  stronger  component  than  would  be  expected  to 
hr*  present  in  the  macroscale  data.  Also,  even  if  the  pack  ice  could  Ik*  considered  a homogeneous 
continuum  at  very  small  scale*-  the  nature  of  the  continuum  constitutive  law  might  couple  with 
variations  in  the  weather  systems  to  give  rapid  variations  in  the  deformation.  This  is  especially 
triii*  since  the  mesoscale  array  is  not  at  the  center  of  the  macroscale  array,  but  at  one  corner  of  it. 
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In  order  to  examine  differences  and  similarities  lietween  the  macroscale  and  mesoseale  de- 
formation, and  thus  test  some  of  the  above  hypotheses,  a comparison  of  the  various  deformation 
t le  series  and  their  spectra  was  made.  The  macroscale  deformation  data  were  supplied  by  Thorn- 
ciiK"  and  are  essentia'iy  the  same  data  as  presented  by  him  (Thorndike  1973).  In  the  processing  of 
the  macroscale  data,  Thorndike  employed  a Kalman  filtering  procedure  to  obtain  this  position  and 
velot  it y of  each  of  t lit*  three  satellite  sites.  The  filter  cutoffs  varied  somewhat  but  generally  all 
frequencies  up  to  about  18-hour  wavelengths  were  passed.  The  strain  rates  and  vorticities  were 
uniquely  determined  since  there  were  only  three  stations. 

Time  series  comparison 

In  order  to  smooth  both  time  series  in  the  same  manner,  all  deformation  rates  were  smoothed 
wit li  a lowpuss  filter  with  a transition  hand  from  21  to  84  hours.  This  smoothing  also  allowed  a 
mesoseale  vorticity  time  series  to  he  constructed  by  including  the  camp  rotation,  a step  that  was 
difficult  without  extensive  smoothing  because  azimuthal  measurements  were  typically  made  only 
once  a day.  To  obtain  the  camp  rotation  for  the  mesoseale  calculation,  a time  series  R (sin<i>) 
longitude  - azimuth  (with  longitude  increasing  in  an  easterly  direction,  see  Nye  (1974)1  with  6 the 
latitude  was  constructed  by  linearly  extrapolating  the  satel'ite  position  and  azimuth  measurements 
reported  hy  Thorndike  et  al.  M972).  Tins  time  series  was  smoothed  by  using  the  same  filter  that 
was  applied  to  the  deformation  data  and  differentiated  to  obtain  a rotation  rate.  This  rotation  rate 
was  added  ti  the  vorticity  calculation  in  tin1  camp  coordinate  system  to  obtain  the  true  vorticity 
as  discussed  in  a previous  section.  For  mesoseale  data  in  the  camp  coordinate  system,  tlie  outer 
array  least  squares  results  were  used.  In  addition,  several  days  of  earlier  data  (taken  before  the 
array  was  complete),  consisting  of  only  4 targets  (1,  2,  3,  8)  and  the  center  point,  wore  used  to 
extend  the  time  series. 

Figure  7 gives  a compaiison  of  mesoscal  and  macroscale  divergence  rates,  vorticities  and 
maximum  shear  rates.  Because  of  a malfunction  of  one  of  the  satellite  navigation  units,  there  is 
a several-day  gap  in  the  macroscale  data  which  was  bridged  by  Thorndike  (1973)  using  a Kalman 
filter.  This  gap  is  denoted  hy  vertical  lines  enclosing  interpolated  macro  data  in  Figure  7.  For 
a quantitative  comparison  of  the  curves,  we  give  in  Table  IV  correlation  coefficients  and  root-mean- 
square  (rms)  values  for  the  various  curves  up  to  the  gap  in  the  macroscale  data.  The  standard 
error  for  the  correlation  coefficients  is  based  upon  a number  of  degrees  of  freedom  equal  to  the  num- 
ber of  points  correlated  times  the  fraction  of  the  spectrum  passed  by  the  filter. 

Figure  7 and  Table  IV  indicate  that  there  are  significant  correlations  between  the  deformation 
time  series  measured  at  different  scales.  Visual  examination  of  the  curves  suggests  that  the 
correlation  is  due  to  the  presence  of  similar  strain  events  over  several-day  periods.  However, 
since  these  events  are  often  of  different  amplitudes  and  occur  at  slightly  different  times,  the  corre- 
lation is  not  complete,  especially  at  higher  temporal  frequencies.  The  results  also  show  that  the 
deformation  rates  have  comparable  amplitudes  with  the  mesoseale  amplitude  generally  being  slight- 
ly larger.  The  comparison  generally  indicates  that  both  meso-  and  macroscale  arrays  are  measuring 
similar  continuum  motions  of  the  ice  pack. 


Table  IV.  Meso-  and  macroscale  root- mean-square  deformation 
rates  and  correlation  coefficients  (from  Fig.  7). 


Divergence  rate  Vorticity Maximum  shear 

Meso  Macro  Meso  Macro  Meso  * Macro 


Root-mean-square  2.84  2.28  4.37  2.45  3 73  2 37 

( 10  ^ hr"1) 


Correlation 

coefficient 


0.47  1 0.20  0.47  1 0.20  0.71  4 0.20  0.71  1 0.20  0.38  * 0 20  0.38  * 0.20 
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Figure  7.  Comparison  ol  mesoscaie  and  macroscale 
divergence  rates,  vorticities,  and  maximum  shear 
rates.  T/ie  dashed  lines  represent  macroscale  data 
and  the  solid  lines  mesoscaie  data. 


Spectral  densities 

Some  of  the  important  differences  between  the  nature  of  the  mesoscaie  deformations  and  that 
of  the  macroscale  deformations  are  illustrated  by  the  spectral  densities  which  we  calculated  using 
the  lap  product  method  as  discussed,  for  example,  by  Rayner  (1971.  p.  94).  Figure  8a  shows  the 
spectra  of  the  mesoscaie  divergence  rate  and  shear  rate  (not  maximum  shear),  and  Figure  8b  shows 
the  spectra  of  the  macroscale  divergence  rate,  shear  rate  and  vorticty.  Because  of  the  inade- 
quacy of  the  camp  rotation  time  series  at  higher  frequencies,  it  is  not  possible  to  construct  a 
mesoscaie  vortie.tv  spectrum.  Also,  because  of  the  Kalman  filter  smoothing  of  the  macroscale  data, 
the  macros,  ale  spectra  are  valid  only  up  .0  about  15-hour  wavelengths.  Since  there  were  M em, 
data  gaps,  the  spectra  did  no,  come  from  the  same  time  periods,  bu,  were  calculated  bom  Julian  day 
8K  to  113  for  the  mesoscaie  data,  and  from  Julian  day  81  ,0  101  for  the  macroscale  data 

Figure  8 shows  tha,  the  macroscale  spectra  generally  contain  less  variance  a,  higher  frequences 
than  the  mesoscaie  spectra.  Such  a result  is  commensurate  with  viewing  the  deformation  as  a co.i- 
tiuiuuu  signal  and  a fluctuation  component  with  the  fluctuation  magnitude  dropping  off  inversely  with 
the  size  of  the  array.  This  follows  because  the  fluctuation  signal,  being  of  a random  nature,  won 
be  expected  ,0  have  a greater  high-frequency  variance  than  the  continuum  signal 

An  interesting  aspect  of  the  mesoscaie  spectra  (especially  the  divergence  rate)  is  the  presence 
of  a significant  spectral  peak  at  about  12-hour  wavelengths.  Whether  such  a peak  is  contained  in 
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the  macroscale  data  cannot  be  ascertained  because  the  smoothing  employed  by  Thorndike  (1973) 
effectively  f.ltered  out  any  such  oscillation.  A possible  explanation  for  this  peak  is  a variation  m 
water  currents  due  to  inertial  oscillations  (Hunkins  1967).  Measurements  or  ocean  eurren  s by 
Newton  and  Coachman  (1973)  during  previous  AIDJEX  pilot  stud.es  indicated  12-hour  cycles  in 
currents  with  the  oscillation  displaying  coherence  over  distances  up  to  20  km.  Different  g 

difteen,  leu  floes  could  couple  wi„,  ,heSc  cum*  » c,ea.c  a differed.,.,  « — ■ 
There  is  also  evidence  of  an  24-hour  cycle  in  the  macroscale  divergence  rate  spectrum  (wine  may 
possibly  be  present  in  the  n.esoscale  data).  Tne  source  of  this  peak  is  not  at  present  understood. 

The  general  fall-off  of  the  spectra  in  Figure  8 is  also  relevant  for  sampling  considerations. 

The  of  the  spool,,  peueralfy  suppes,,  ,.,a,  Samp„ne  intervals  up  to  .0  hours  (mb  acc«e 
measurements)  would  yield  low  frequency  information  without  intolerable  aliasing.  A more  dir 
test  of  this  can  be  made  by  sampling  the  data  at  larger  intervals  before  smoothing  and  comparing 
results  to  data  smoothed  before  resampling.  Such  comparisons  have  been  made  for  the  mesoscale 
data  (Hibler  et  al.  1973b)  and  support  the  conclusion  that  accurate  samples  every  8 hours  are  ade- 
quate for  resolving  the  low  frequency  components  of  the  time  series  required  for  comparison  with 
synoptic  Meteorological  variations  which  generally  occur  over  a time  scale  of  several  days  (Monin 

1972). 


Nature  of  the  Ice  Pack  Rotation 


Examination  of  the  mesoscale  vorticity  indicates  that  it  is  similar  to  the  camp  ro  at  ion.  This 
can  be  seen  from  Figure  9,  which  shows  the  camp  rotation  rate  and  the  mesoscale  vorticity.  This 
similarity  means  that  to  a large  degree  the  whole  mesoscale  region  is  rotating  as  an  entity.  Inves  i 
gation  of  the  macroscale  deformation  data  indicates  that  such  a "solid  rotation  is  also  Partially 
o« " ring  for  The  larger  macroscale  region,  a,  least  at  low  temporal  frequencies.  This  may  be  seen 
from  Figure  10  in  which  the  east-west  shear  rate  «V‘1x  and  vorticity  are  plotted  together.  As  can 
be  seen,  both  curves  are  similar,  indicating  that  *x/*y  and  are  close  to  being  equal 


IS 


DIFFERENTIAL  SEA  ICF  DRIFT 


. o oi2«  1 *0*§M»*  Wa1* 

94  " «8  * * 5*  ‘ i® 

T.«W*,  dOf» 


Figure  9.  Camp  rotation  rate  and  me.sosca/e  vorticity. 
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Figure  10.  Comparison  of  east-west  shear  rate  (<lvy 
Ox)  and  vorticity. 

magnitude  and  opposite  in  sign,  a condition  which  holds  identically  for  solid  rotation.  Another 
indication  of  this  behavior  on  the  macroscale  is  the  relative  smallness  at  low  frequencies  of  the 
macroscale  shear  rate  spectral  density  in  the  x-y  coordinate  system  compared  with  the  macroscale 
vorticity  spectral  density  as  shown  in  Figure  8b. 

These  observations  suggest  that  the  most  pronounced  mode  of  differential  deformation  of  the 
ice  pack  is  a relatively  cohesive  rotation,  at  least  at  low  frequencies,  of  the  pack.  The  direction 
of  the  rotation  is,  from  the  curves  in  Figure  4,  generally  clockwise  in  agreement  with  the  motion 
of  the  Pacific  Gyre.  Other  deformation  rates  at  low  frequencies  appear  to  be  somewhat  smaller 
than  the  rotation  rate.  Such  a cohesive  rotation  is  also  borne  out  by  shear  zone  deformation  studies 
(Crowder  et  al.  1974,  Hibler  et  al.  1974),  suggesting  that  the  pack  is  rotating  as  a relatively  tightly 
bound  continuum  with  slippage  at  the  boundaries. 
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(see  °f  *7?  ,S  “ n('"at,V(>  ('Orr0li,tl0,1  vort icily  and  the  divergence  r ,,«■ 

set  fig.  7).  lie  correlation  coefficient  for  the  n,esoscale  data  (omitting  the  earlier  d ..a  • 

b(  lore  the  njesosciile  array  was  complete)  is  0 1rH  fi  17 
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Conclusions 

We  believe  that  this  study  has  shown  that  strains  measured  on  1 se'ita  in  i-, 
serve  as  a valid  measure  of  the  continuum  motion  of  sea  ice.  When  strains  are  me-isurerT^- * 1 M 
smaller  scale  the  continuum  motion  is  obscured  by  the  random  bumping  and  viewing  of  floes’and 
by  he  finue  size  o the  interacting  floes  relative  to  the  size  of  the  strain  army,  ll  se  s^  .l 

ful  inrr  f ' CVer’  in,eres,inK  1,1  th(‘u  own  r*K'"-  Ota  r, -suits  also  indicate  that  me  ming- 
indhS, CU"  b“  fl0ni  ^ **  tnea,„r,d  rota, ‘ton  of 

These  results  suggest  that  it  should  be  possible  to  carry  ou,  highly  useful  experiments  on  the 
mm-roscale  (continuum,  behavior  of  pack  ice  by  utilizing  one  manned  daftin^ZiT C ” ih 
• la.  transponders  located  a,  the  remote  strain  points.  S.rnn  arrays  with  effective  diameters  of 
1 to  nO  kin  should  easily  oc  possible  using  currently  available  techniques. 

F.naHy.  it  should  he  emphasized  that  our  results  were  obtained  from  a pac  k ice  field  with  a 

r^wnLtX0n^a^  U T**'  b,rt  Unkn°WUl  "*  ,h‘Pkw*«  ciis.rilin, ion.  Whether  similar 
W0Uld  b(  ob,,nnC(1  lor  we  conditions,  is  presently  a moot  point. 
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PART  n.  COMPARISON  OF  MESOSCALE  STRAIN  MEASUREMENTS 
WITH  LINEAR  DRIFT  THEORY  PREDICTIONS 


by 

W.D.  Hibler  111 


Introduction 


One  of  the  more  important  uses  of  mesoscale  strain  measurements  is  the  comparison  of  the 
differential  drift,  i.e.,  strain  results,  with  sea  iee  drift  theories,  both  to  test  the  theories  and  to 
determine  eertain  unknown  parameters.  These  comparisons  provide  a more  critical  measure  of 
certain  constitutive  law  parameters  than  do  comparisons  involving  the  drift  of  only  a single  point. 
To  make  such  a comparison  in  this  report,  a linearized  drift  theory  similar  to  that  used  by  Egorov 
(1970,  1971),  Rothroek  (1972),  and  Witting  (1972)  will  be  used.  Such  theories,  although  not  as 
exact  as  other  calculations  (Campbell  1965,  Campbell  and  Rasmussen  1973),  do  suffice  for  obtain- 
ing quantitative  estimates  of  the  dominant  drift  effects. 

Recent  calculations  using  line  r drift  theories  have  generally  been  of  two  types.  In  the  first, 
as  earried  out  by  Rothroek  (1972)  and  Witting  (1972)  independently,  the  average  yearly  circulation 
of  the  arctic  iee  cover  has  been  calculated  assuming  that  the  iee  is  incompressible.  Such  calcula- 
tions, although  of  considerable  interest,  are  not  directly  comparable  to  mesoscale  strain  measure- 
ments. In  the  seeond,  as  carried  out  by  Egorov,  an  approximate  infinite  boundary  solution  is  ob- 
tained to  a linear  drift  theory  that  uses  a shear  viseosity  to  explain  the  rheological  behavior  or  the 
iee  and  neglects  gradient  eurrent  effects,  i.e,,  neglects  geostrophic  ocean  flow  below  the  surface 
boundary  layer.  The  second  approach  yields  results  more  directly  comparable  to  strain  measure- 
ments. 

For  the  comparison  made  here,  a linear  drift  equation  similar  to  that  used  by  Egorov  will  be 
used.  The  rheological  behavior  of  the  iee  will  be  taken  into  aceou.it  by  using  the  constitutive  law 
proposed  by  Glen  (1970),  whieh  includes  a bulk  viseosity  as  well  as  a shear  viscosity.  In 
addition,  calculations  will  be  made  using  a more  general  linear  viseoelastie  constitutive  law  that 
allows  for  memory  effects  and  includes  a generalized  Hooke’s  law  as  well  as  the  Glen  law  as 
special  cases.  As  in  Egorov’s  work,  gradient  current  effects  will  be  neglected. 

However,  unlike  Egorov,  we  will  formulate  the  infinite  boundary  solution  without  approxima- 
tion in  terms  of  a linear  response  function.  The  resulting  real  space  solution  consists  of  a 
straightforward  integral  operator  which  may  be  applied  to  the  pressure  field  to  obtain  expected 
differential  ice  drift.  This  response  function  form  of  the  solution  is  useful  because  it  clearly 
illustrates  the  differences  of  ice  drift  behavior  expected  in  winter  as  opposed  to  that  expeeted  in 
summer.  In  addition,  the  solution  allows  a rapid  determination  of  the  scales  of  variation  in  the 
atmospheric  pressure  field  that  are  important  for  given  bulk  and  shear  viscosity  values.  Further- 
more, the  comparison  of  strain  measurements  with  pressure  data  allows  the  bulk  and  shear  viseosity 
parameters  to  be  estimated 


List  of  Symbols 


w ice  vorticity 


sec 


A iee  divergence  rate 


see 


differential  sea  ice  drift 
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ice  mass  per  unit  area 

kg  m 2 

11 

ice  velocity 

m sec-1 

V 

wind  velocity 

m sec  1 

t 

Coriolis  vector 

see-1 

t 

C oriolis  parameter  equal  to  magnitude  of  the  Coriolis  vector  times 
the  sine  of  the  latitude 

sec  1 

\ 

mf 

kg  see"1  nf~ 

F 

force  due  to  internal  ice  stress 

N m ’ 

'-w 

water  stress  on  ice 

N m ~ 

T 

a 

an  stress  on  ice 

N in"' 

«/> 

Hknian  angle  in  air 

dimensionless 

0 

Ekman  angle  m water 

dimensionless 

B 

wind  stress  constant  f>(IK,/2)'‘ 

kg  sec  1 m ~ 

D 

water  stress  constant  n (IK  /2)' 

kg  sec  1 m ~ 

/> 

air  density 

kg  nf 3 

/’w 

water  density 

kg  m 

K 

a 

eddy  viscosity  of  air 

;>  1 
ill  sec 

K 

w 

eddy  viscosity  of  water 

in  sec  1 

11 

it 

x component  of  geostrophic  wind 

m sec"1 

V 

y component  of  geostrophic  wind 

m sec ' 1 

P 

atmospheric  piessurt* 

mli 

P 

spatial  average  ot  the  atmospheric  pressure 

mb 

shear  viscosity  of  ice 

kg  sec" 

4 

bulk  viscosity  of  ice 

kg  sec" 

H 

divergence  rate  response  function 

dimensionless 

G 

vorticity  response  function 

dimensionless 

Linear  Drift  Equations 

h ol lowing  Egorov  s (1970  1971)  example,  we  consider  a steady  state  equilibrium  drift  eqna- 
ton  for  the  case  where  the  gradient  current  term  varies  so  slowly  in  space  and  time  that  it  may  be 
u g e<  ec . y neglecting  this  term,  the  sea  ice  is  effectively  considered  to  be  moving  across  a 

stagnant  ocean  [see.  for  example.  Rotluoek  (1972)1.  In  this  case  the  equilibrium  equation  takes 
the  form 


-mf  y u r rw  1 ra  t F 


0 


(1) 


where  u is  the  ice  vekx  ity,  / tne  Coriolis  vector,  m the  ice  mass  per  unit  area.  F the  force  due  to 
internal  ice  stress,  and  rw  and  ra  the  water  and  air  stresses  on  ice,  respectively.  The  components 
of  water  and  air  stresses  are  given  by  a simple  Ekman  layer  theory 
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Ta  6 (cos  £>U  sin(.6F,) 

X " K 

„ (2) 
ra  B (*  cos  <j>  V * sin  6 U ) 

V h K 

rw  D (-  cos  0 ii  . sm  0 u ) 

x * V 

\ - D (sin  0 it  t cos  0 u ) 

y * V 

where  6 mid  0 are  the  Kkman  angles  in  the  air  and  water,  respectively.  The  parameters  B and  D 
are  proportionality  constants  related  to  the  turbulence  coefficients  for  the  atmosphere  and  the 
ocean.  For  the  classical  Ekman  layer  solution  (Sutton  19r>3,  p.  71),  B and  D are  given  by  p\(ll<:2) 
where  p and  K are  the  air  density  and  eddv  viscosity  respectively  of  either  an  or  water.  U and 
itr*’  'lie  geostrophic  wind  components  given  bv  w 


U 

K 


1 AP 
I'l  By 


V 

K 


1 r IP 
pi  Ax 


(da) 


(lb) 


where  p is  the  an  density.  P is  the  atmospheric  pressure  and  I is  the  magnitude  of  the  Coriolis 
vector.  In  eq  2 it  is  implicitly  assumed  that  the  ice  velocity  is  small  compared  with  the  wind 
velocity  and  may  be  neglected.  For  F , the  force  due  to  internal  ice  stress,  the  constitutive  law 
proposed  by  Glen  (1970)  is  used 

f v^ut<r(r  u)  (5) 

where  £ and  rj  are  bulk  and  shear  viscosity  constants  that  can  vary  with  ice  compactness  and  thete- 
fore  season.  Calculations  using  a somewhat  more  general  constitutive  law  are  discussed  later. 


Ice  Drift  Solutions 

In  this  case  the  interest  is  primarily  in  the  solutions  of  the  linear  drift  equations  for  the  ice 
divergence  rate  A I \ (Au%  Ax)  < (Au  /Ay)\  and  ice  vorticity  w|w  l/2[(Au  /Ax)  - (Au  /By)]|,  By 
taking  the  divergence  and  curl  of  eq  1,  two  linear  equations  for  A and  w are  obtained  * 

[(17  t Or"’  - D cosO]A  4 (A  * D sinl?)2w  B Sin<^  P (6) 


-(A  » D sin0)A  , (r/^~  - D costf)2 w — cos_  «!-  p 

pl 


whete  A ml  These  equations  represent  a linear  system  with  the  input  being  the  pressure  field  P 
and  the  output  being  A and  w.  Such  systems  (Jenkins  and  Watts  1968)  may  be  described  by  response 
functions  in  wave  number  space  H ,(k).  //,(*).  so  that  A (k)  H,(k)  P(k ) and  w(k)  HJk)  P(k)  where 

we  denote  wave  number  space  functions  with  a tilda  and  k k . The  response  functions  mav  be 
obtained  straightforwardly  by  Fourier  transforming  eq  6 and  7,  yielding  wave  number  space  equa- 
tions 
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\( k)  B-  — ll  H(k)  1 

pl  'I  » <= 


(8) 


w(k)  — — - 1 1 G(k)l 
pl  2 »/ 


(9) 


where 


1 II  (k)  — ° c'os(!,'  S‘»  $ ~ cos  ^ (A  1 P sin  ^ 

IX2  * D2  . 2D  sin  CM  » (t/  t <)  r^k4  4-  D cosO  (2 7 4 Ok2) 

k2il(r,  . Ok'  * D cosOl  cos</>  - sin £>  (A  + D sinWlq 
1 Gtk)  ■ . 

U2  . D2  ♦ 2D  sin  f>A  4 (7  4 £)  17k4  4 D cos  0 (2r/  . £)k"l 


(10) 


(ID 


By  the  convolution  theorem  these  equations  yield  simple  integral  equations  in  real  space.  For 
example,  for  \(x) 


\(x)  - 


B 


pl(rj  4 £) 
where  (using  polar  coordinates) 


IP(x)  P'(x)| 


(12) 


P'(x)  4 j IH r - r )f  Pit'.  0)  rdOdv 


(13) 


and 


//( r)  i-  f !Dk)kJQ  (kr)dk.  (14) 

2n  0 

and  J()  is  the  zero  th  order  Bessel  function.  It  should  be  noted  that  t^ese  equations  only  apply 
exactly  to  an  ice  cover  and  pressure  field  of  infinite  extent  and,  in  fact,  represent  solutions  using 
boundary  conditions  P(x),  \<x),  w(x)  finite  at  x,  y • ♦ However,  in  practice  they  may  be  applied 
to  a finite  case  with  the  necessary  extent  of  the  ice  cover  and  pressure  field  determined  by  the 
spatial  extent  of  a finite  filter  W(r)  that  approximates  the  wave  number  response  H(k)  in  eq  10. 

The  wave  number  space  form  of  the  response  functions  IKk)  and  G(k)  contains  considerable 
information.  In  the  case  of  the  divergence  rate  for  q,  £ large,  1 - H(k)  is  generally  positive  for 
large  k and  negative  or  zero  for  small  k.  Thus,  the  divergence  rate  is  essentially  the  result  of  a 
high  pass  filtering  operation  on  the  pressure  field  with  the  high  wave  numbers  contributing  posi- 
tively to  the  divergence  rate  and  the  low  wave  numbers  negatively  with  a smaller  amplitude.  For 
the  vorticity,  the  response  function  is  also  a high  pass  filter,  but  there  is  no  change  in  the  sign 
of  the  contribution  from  different  wave  number  components  of  the  pressure  field.  The  wave  number 
of  the  filter  cutoff  decreases  as  q and  £ increase.  Thus,  different  types  of  behavior  are  expected 
for  7,  £ small  as  compared  with  7.  £ large.  These  different  types  of  behavior  may  be  characterized 
by  examining  the  two  limiting  cases  t),  £ « 0 and  ij,  £ * 00  • 
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Limiting  cases 

For  tin*  first  case,  noting  that  |1  - H(k)|/(q  4 Oi2  and  ll  - G(k)l  qk2  are  finite  as  q,  £ - 0. 
we  have 


Lini 

V.£  *0 


A (1) 


-B  2 . lOlcosW  sin</>  - sin#  cos<£)  - A cos<Al 

— ^ P(x) — 

^ A2  4 D2  4 20A  sin  0 


(15) 


Lim 
Ti.C  -0 


wtx) 


— r2  p>*> 

2pf 


iDtcosO  cos  d>  • sin0  sin<£)  4 A sincil 
A 2 . D2  * 2DA  sin  0 


(16) 


Thus,  we  see  the  well  known  result  that  for  equal  Ekman  angles  and  small  q,  £ the  ice  would  be 
expected  to  diverge  in  a low  (C2P  * 0)  and  converge  in  a high,  whereas  the  vorticity  would  be 
expected  to  be  positive  (counterclockwise)  in  a low  and  vice  versa  in  a high. 

For  the  second  limiting  case,  we  note  that,  for  q,  £ very  large.  H(k)  and  G(k)  pass  only  the 
verv  long  spatial  wavelengths  with  the  pass  band  frequency  cutoff  scaling  as  1/q  and/or  l/£ . 
Consequently,  the  real  space  response  functions  H(x>  and  G(x)  approach  constants  (with  integrated 
areas  of  unity)  for  very  large  q.  £.  Therefore.  /W(x  x')P(x')  dx'  and  I G<x  x’)P(x')  dx‘  approach 
the  average  pressure  for  large  q and  £:  as  a result  the  large  q,  £ limiting  equations  are 


Lim 

A (x) 

IPlx)  - P I sin</> 

(17) 

!■£  •» 

(q  . t,  )f 

Lim 

W(x)  - 

— 1 P(x)  Pleostl 

(18) 

q.  £ •» 

2 qf>( 

when-  P is  the  mean  atmospheric  pressure  oyer  the  infinite  x.  y space  which  would  be  approximately 
constant  in  time.  For  1/.  £ large  but  finite.  P would  he  replaced  by  the  very  low  wave  number 
components  of  the  pressure  field  which  would  be  expected  to  he  reasonably  constant  in  time  if  the 
cutoff  wavelength  were  longer  than  the  synoptic  variation  scale  of  the  pressure  field. 

As  can  be  seen  from  eq  17  and  18.  in  the  large  17.  c limiting  case  the  divergence  rate  and 
vorticity  are  proportional  to  the  local  pressure  deviation  from  the  overall  mean  pressure  with  a 
low  pressure  indicating  a convergence  and  a positive  vorticity.  Note  that  there  is  no  dependence 
on  the  water  stress  in  this  limiting  case.  In  fact  the  large  q.  £ case  is  equivalent  to  neglecting 
all  stresses  except  the  internal  ice  stress  and  wind  stress.  An  alternative  derivation,  for  example, 
would  be  to  delete  the  water  stress  and  Coriolis  terms  from  eq  6 and  7 and  solve  a boundary  value 
problem  with  (A  - P)  and  (w  • P)  finite  at  x.  y • ♦ * . It  is  also  important  to  note  that  the  large  _ 

1/,  £ solution  includes  lateral  transfer  of  stress  through  the  pack  up  to  infinite  distances  via  the  P 
term.  However,  this  term  becomes  only  a constant  because  the  lateral  stress  averages  out  and 
thus  A and  w follow  the  local  pressure. 

To  the  extent  that  q.  £ may  he  considered  very  large  in  the  winter  ami  small  in  the  summer,  the 
two  limiting  cases  suggest  that  sea  ice  (far  from  coastal  boundaries)  would  be  expected  to  converge 
in  a low  in  winter  and  diverge  in  a low  in  summer  with  vorticity  always  positive  in  a low.  Such 
predicted  behavior  agrees  with  earlier  mesoscale  strain  measurements  (Hibler  et  al.  1973).  with  the 
more  extensive  results  reported  in  this  report,  and  with  Russian  observations  (Volkov  et  al.  1971). 

It  is  also  what  one  would  expect  intuitively,  namely,  that  in  winter  the  ice  is  tightly  held  and  can- 
not move  rapidly,  so  that  the  water  and  Coriolis  forces  would  be  expected  to  he  smaller  than  in 
tlie  summer. 
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figure  1.  Wave  number  space  response  functions  for  the  divergence  rate  and  vorlicity  of  the  ice 
puck  for  different  values  of  and  c,.  The  response  functions  operate  on  the  atmospheric  pressure 

field. 


F inally,  we  note  that  it  wo  used  a series  solution  for  the  drift  as  given  by  Egorov  (1970)  it 
would  be  impossible  to  draw  the  above  conclusions,  because  the  series  diverges  for  frequencies 
higher  than  the  high  pass  cutoff  frequency  in  1 - Hik). 

General  case  - wavelength  dependence 

Clearly,  it  is  critical  what  the  "cutoff  wavelengths”  for  1 - H(k)  and  1 - G(k)  are.  To  illus- 
trate typical  forms  ot  lUk)  and  G(Jc),  we  use  the  following  numerical  values: 


f 

1.46 

■ 10 

4 sec' 1 

m 

3.0  • 

10s 

kg  m 3 

0,6 

30 

l> 

1.3  kg  m 

3 

“a 

1.5  ■ 

10 1 

m"  sec  1 

2.0  • 

10  ‘ 

J »>  1 

“ nr  sec' 

B 

4.3  • 

10  kg  sec  1 

D 

1.18 

1 1 

kg  sec  m . 

The  results  for  different  values  of  7 am  C illustrated  in  Figure  1.  For  the  H(k)  curve  the  key 
wavelength  is  the  transition  from  positive  to  negative  response.  The  results  generally  indicate 


MESOSCALE  STRAIN  MEASUREMENTS  AND  LINEAR  DRIFT  THEORY  PREDICTIONS  27 


that  for  7 , £ 10l“>  kg  sec-1,  the  "high  wave  number"  pressure  variations  are  those  with  wave- 

lengths shorter  than  11000  kin.  Since  the  synoptic  scale  for  the  Arctic  is  of  the  order  of  500  km 
(Egorov  1971),  the  large  7,  C limiting  case  is  expected  to  lx*  reasonably  valid  for  7.  £ 10*^  kg 

sec  ' . For  7,  £ lo'1  kg  sec  ',  on  the  other  hand,  the  high  wave  number  cutoff  is  such  that  one 
might  expect  the  small  7,  £ case  to  be  more  applicable,  This  wavelength  dependence  explains 
why  most  calculations  by  Campbell  (1965)  and  Campbell  and  Rasmussen  (1973)  have  indicated 
diverging  ice  in  a low  press  ire  region.  This  is  especially  true  of  the  yearly  average  drift  where 
the  mean  yearly  pressure  fie  d contains  few  high  wave  number  spatial  variations. 

Figure  1 also  illustrates  the  importance,  of  the  scale  over  whir'll  the  pressure  field  is  consider- 
ed. For  example,  consider  7 and  Rvalues  such  that  the  positive  to  negative  transition  111  H(k)  is 
at  about  1000  km.  Then,  considering  a high  pressure  system  varying  slowly  in  space  with  few  high 
frequency  components  with  wave  numbers  greater  than  1000  km-1  the  ice  would  be  expected  to 
converge.  On  the  other  hand,  for  a high  pressure  system  varying  very  rapidly  in  space  with  signifi- 
cant variance  at  wavt  numbers  greater  than  1000  knr  *,  the  ice  would  b-*  expected  to  diverge  where 
the  pressure  was  high.  Clearly,  it  is  very  important  to  define  the  spatial  scale  used  when  speaking 
of  diverging  or  converging  ice. 


Comparison  of  Theory  with  Mesoscale  Measurements 

To  determine  how  well  the  limiting  forms  of  the  predicted  A and  w values  for  large  7 and  £ and 
infinite  boundaries  compare  with  mesoscale  observations  reported  in  Hibler  et  al.  (1974),  the  local 
pressure  at  the  main  A1DJEX  1972  Camp  (located  at  roughly  75  N 148  W),  the  measured  ice 
divergence  rate,  and  the  measured  ice  vort icily  were  compared.  The  resulting  time  series  are 
illustrated  in  Figure  2 with  the  dashed  portions  of  the  deformation  rates  representing  data  taken 
while  the  mesoscale  array  was  only  partially  deployed.  Calculation  of  the  strain  rate  and  vorticity 
time  series  is  described  in  some  detail  in  Part  1 of  this  report.  In  addition  to  these  three  time  series, 
Figure  2 shows  the  calculated  divergence  of  the  w»nd  velocity  field  and  the  fluctuations  of  the 
local  atmospheric  pressure  from  the  average  pressure  over  an  approximately  600-km-diameter  region. 
The  average  pressure  P was  estimated  by  taking  the  average  of  the  camp  pressure,  four  remote  data 
buoy  pressures  located  around  the  camp  about  300  km  away,  and  the  Point  Barrow  pressure.  For 
calculation  of  divergence  of  the  wind  velocity  field,  local  wind  speed  and  direction  measurements 
at  each  of  the  three  manned  stations  were  used.  The  distances  and  relative  angles  between  the  sta- 
tions were  taken  as  constant  and  estimated  from  position  data  for  19  March  as  reported  by  Thorn- 
dike et  al.  (1972).  The  basic  computational  equations  are  similar  to  those  use"1  in  the  strain  cal- 
culations in  Part  I of  this  report.  All  of  the  time  series  shown  in  Figure  2 were  smoothed  with  the 
same  low  pass  filter  having  a transition  band  from  0-3/80  cycles  hr  1 (Hibler  19~2). 

Correlation  coefficients  were  calculated  between  all  five  of  the  time  series  (excluding  the 
dashed  portions  of  the  divergence  rate  and  vorticity)  with  the  results  listed  in  matrix  form  in  Table 
1.  The  standard  error  is  based  upon  a number  of  degrees  of  freedom  equal  to  the  number  of  points 
correlated  times  the  fraction  of  the  spectrum  passed  by  the  filter. 

As  can  be  seen  from  Figure  2 and  from  the  table  of  correlation  coefficients,  there  is  a positive 
correlation  between  the  local  pressure  and  the  divergence  rate  and  a negative  correlation  between 
pressure  and  vorticity  as  predicted  by  the  large  viscosity  limiting  case  of  the  linear  drift  theory. 

The  results  also  indicate  that  the  spatial  pressure  fluctuation  time  series  is  quite  similar  to  the 
camp  pressure  time  series  and  has  a similar  correlation  to  the  divergence  rate  and  vorticity.  This 
generally  indicates  that  the  pressure  field  has  cor  siderable  variance  at  wavelengths  shorter  than 
600  km  and  justifies  to  a limited  extent  the  use  of  the  infinite  boundary  solution  for  comparison. 
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were  smoothed  with  a low  pass  tiller  having  a transition  hand  trom  0 to  3/80  cycles  hr-1. 


Table  I.  Correlation  coefficient  matrix  between  time  series. 

u ice  velocity,  v wind  velocity  (standard  error  0.22). 


V » 

V x lt/2 

P - P 

P 

v-  y 

V-  « 

1 

-0.44 

0.41 

0.34 

0.38 

C v » '2 

-0.44 

l 

0.54 

-0.46 

-0.53 

p -P 

0.41 

-0.54 

1 

0.75 

0.39 

p 

0.34 

-0.46 

0-75 

1 

0.21 

r-  y 

0.38 

-0.53 

0.39 

0.31 

1 
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Frequency  (cycles  hr*') 

Figure  3.  Spectra  ol  atmospheric  pressure  and  me s<>- 
scale  ice  divergence  rale  at  the  main  AIDJEX  1H72 
camp. 

Thi*  correlation  between  the  vorticity  and  the  divergence  rate  shows  the  expected  negative 
value,  with  the  magnitude  of  the  vorticity  being  generally  larger  than  that  of  the  divergence  rate. 

In  partic  ilar.  the  ratio  of  the  vorticity  variance  to  the  divergence  rate  variance  is  :}.8.  The 
correlations  between  the  wind  divergence  and  the  local  pressure  and  pressure  fluctuations  are 
also  positive,  indicating  the  expected  wind  convergence  in  a low  pressure  region  and  vice  versa 
in  a high  pressure  region. 

With  respect  to  correlation  at  higher  temporal  frequencies,  then*  are  indications  that  there  is 
little  linear  correlation  between  the  pressure  and  the  divergence  rate  at  periods  shorter  than  24 
hours.  This  is  reasonable,  since  the  atmospheric  pressure  variation  is  veiy  nearly  band  limited. 
This  is  illustrated  by  the  spectia  of  the  pressure  divergence  rate  («n)  time  series  given  in  Figure  2. 
Clearly  the  pressure  time  series  has  comparatively  little  variance  at  periods  shorter  than  '.'4  hours, 
which  agrees  with  typical  expected  synoptic  variation  scales  (Monin  1972.  p.  9).  The  divergence 
rate  spectrum,  on  tlie  other  hand,  is  relatively  flat,  although  it  does  fall  off  bv  a factor  of  about  2 at 
24-hour  periods.  These  curves  indicate  that,  although  the  meteorological  driving  forces  cm  the  ice 
are  relatively  smooth,  tin*  response  of  the  ice  is  more  complex  and  erratic  hi  time  probably  due  to 
random  bumping  of  flows  and  opening  and  closing  of  leads. 

Estimates  of  coMtitsttve  law  para  alters  rj  and  £ 

Assuming  that  the  large  viscosity  limiting  case  of  the  linear  drift  theory  is  applicable  to  out 
observed  mesoscale  deformation,  we  may  estimate  the  viscosity  parameters  </  and  c . To  do  this  we 
utilized  the  slopes  of  the  regression  lines  of  A and  w upon  P (using  the  curves  m Figure  2).  and 
also,  for  comparison,  the  regression  lines  of  A and  w upon  P - P.  F.quating  these  regression  line 
slopes  to  the*  predicted  slopes  in  eq  17  and  18  and  inserting  numerical  values  fix  B.  </.  and  I (as 
previously  listed),  we  obtained  the  results  shown  in  Table  II. 
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Table  D.  Estimate*  of  r,  tad  £ . 

Regression  line  used q ^ t £ 

\ and  w u|>on  P (1.70  i 0.73)  • lO*2  kg  sec"1  (5.3  ♦ 3.34)  • I012  kg  sec-1 

A ami  w upon  P - P (0.66  f 0.23)  • IQ12  kg  sec-1  (1.94  i 0.95)  • 1012  kf  sec  1 


According  to  the  estimates  in  Table  II.  r/  and  £ are  of  the  order  of  1012  kg  sec  1 with  £.  the 
bulk  viscosity,  being  somewhat  larger  than  i),  the  shear  viscosity.  This  result  agrees  with  intuitive 
expectations,  since  one  would  expect  the  ice  to  offer  greater  resistance  to  pure  compression  than 
to  pure  shear. 

Table  11  also  shows  that  the  estimates  of  '/  and  £ from  different  regression  lines  are  quite 
similar.  This  illustrates  that  much  of  the  observed  correlation  between  the  pressure  field  and  the 
differential  ice  motion  is  due  to  higher  wave  number  (spatial)  variations  in  the  pressure  field.  Such 
behavior  suggests  an  explanation  for  why  the  infinite  boundary  solution  works  reasonably  well 
since  high  wave  number  variations  may  be  extracted  by  a real  space  response  function  that  is  well 
limited  in  space.  Response  functions  to  extract  lower  wave  number  variations,  on  the  other  hand, 
extend  much  further  spatially  and  boundary  effects  would  consequently  be  expected  to  be  more 
critical  at  lower  wave  numbers. 

Note  also  that  reducing  B,  the  wind  stress  constant,  would  reduce  q and  £,  and  cause  the 
response  function  to  be  more  limited  in  space.  Keeping  this  factor  in  mind,  it  is  likely  that  our 
estimates  of  q and  £ are  only  approximate  and.  in  general,  our  deformation  results  could  be  com- 
mensurate with  rj  and  £ values  varying  anywhere  from  10M  to  1012  kg  sec-1. 

Regarding  the  validity  of  the  large  viscosity  limiting  case,  by  referring  back  to  Figure  1 we 
see  that  for  values  of  q,  £ 1012  to  1013  kg  sec'1  the  wave  number  cutoff  of  the  response  functions 

lies  in  the  range  1000  to  3000  km.  These  wavelengths  are  commensurate  with  or  larger  than  ex- 
pected synoptic  variations  in  the  pressure  field,  so  that  the  use  of  the  large  q.  £ limiting  case 
appears  to  be  justified  for  the  data  analyzed  in  this  paper.  However,  for  the  smaller  wavelengths  in 
this  range,  the  expected  correlation  would  be  primarily  between  the  high  wave  number  pressure  com- 
ponents such  as  those  estimated  using  P - P.  The  large  q,  £ limiting  case  may  also  be  justified 
for  other  boundary  conditions  (see  App.). 

It  is  interesting  to  note  that  compressive  stresses  predicted  by  our  estimated  values  of  q and 
£ are  reasonable  in  terms  of  stresses  predicted  by  Parmerter  and  Coon  (1972).  For  example,  maxi- 
mum values  of  A are  of  the  order  of  0.0004  hr"1,  which  yields  for  £ 1012  g sec"1  a compressive 

stress  of  1.1  • 10->  N m 1 where  we  have  used  the  Glen  constitutive  law.  This  is  close  to  the  0.1 
to  0.4  • 10'*  N m needed  to  cause  ridging  in  2-m  ice  by  the  bending-failure  Parmerter  and  Coon 
ridge  model.  It  is  also  similar  to  the  maximum  pressure  difference  of  2 ■ 105  N m_1  obtained  by 
Rothrock  (1972)  assuming  the  ice  is  incompressible. 

Also  of  some  interest  is  the  numerical  comparison  of  the  calculated  wind  divergence  rate  with 
that  estimated  from  the  curvature  of  the  pressure  field.  To  do  this  we  take  as  an  estimate  of  the 
Lapacian  - (4/5)  [(P  - P)/a2l  where  a • 300  km.  Using  the  regression  line  of  V-  y upon  P - P, 
we  find  an  observed  relation  which  yields  an  Ekman  angle  of  44°.  Certainly  these  comparisons  are 
only  approximate,  but  they  do  indicate  that  the  wind  divergence  estimated  from  the  pressure  field, 
using  the  geostrophic  approximation  and  a constant  Ekman  angle,  is  the  same  order  of  magnitude  as 
the  calculated  wind  divergence  rate. 
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A More  General  Linear  Constitutive  Law 

Tlit>  previous  sections  generally  indicate  that  most  of  the  dominant  aspec  ts  of  the  observed 
mesoscale  drift  behavior  may  be  explained  using  a simple  viscous  iee  rheology.  It  is  useful  to 
see  if  better  agreement  can  be  obtained  using  a linear  viscoelastic  law  including  me  mory  effects. 
Io  do  this  we  will  carry  out  calculations  using  a more  general  constitutive  law  which  allows  the 
“viscosities"  (bulk  and  shear ) to  vary  with  frequency  and  which  can  include  both  elastic  and 
viscous  behavior.  One  sucli  law  t hat  is  computationally  similar  to  the  Glen  law  is  given  by 


F(t ) 


t )yly  * u(r)idr 


(19) 


where  t denotes  time. 

Taking  the  temporal  Fourier  transform  of  ttiis  equation,  we  obtain  (for  convenience  we  simply 
replace  t by  o>  to  denote  temporal  transforms) 

F(cj)  »/(</i)V"  u (o*)  . c (o)  vl  y • u ((,>)|  (20) 

where  t/(o j)  and  £(oj)  ale  analytic  m the  upper  half  plane  to  guarantee  causality.  Two  particular 
limiting  cases  of  this  law  are 

1)  Glen  viscous  law 


7/(t)  r/S(t) 

£(0  <S(t) 


(31) 


or,  in  frequency  space 
77(01)  77 

CM  C 


(22) 


where  77  and  £ are  constant  viscosities. 
2)  Generalized  Hooke's  law 


where 


77(1) 

7/  0(0 

on 

£0(0 

(HO 

1 t go 
0 t < 0 

(23) 


(24) 


or,  in  frequency  space 
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,l(oj)  ’I 

Oi  + IO 

(25) 

1 ...  £ 

(o  ♦ id 


with  S infinitesimal. 


Note  that,  for  the  Hooke's  law  case,  r/lw)  and  £(u>)  both  decrease  with  decreasing  frequency  and 
have  a phase  shift.  The  phase  shift  is  the  key  indicator  of  elastic  behavior. 

Drift  calculations  using  generalized  law 

It  is  clear  that  by  using  temporal  transforms  of  all  quantities  the  same  formalism  used  for  the 
simple  viscous  calculation  may  be  used  for  the  more  generalized  calculations.  In  particular,  eq  6- 
11  may  be  formally  extended  to  include  the  generalized  constitutive  law  by  replacing  all  quantities 
with  temporal  Fourier  transforms;  for  example,  A(x)  . \(x,  o>);  A(k)  . A (It.  u>).  The  arguments 
about  limiting  cases  also  are  similar  except  that  the  magnitudes  ot  q(u>)  and  £(o>)  are  now  the 
determining  factors.  In  particular,  for  the  large  i»/(ia)  . !£(«)!  limiting  case,  we  have  by  analogy 
to  eq  17  and  18  the  results: 


Lim  A(x, (.»)  BlP(x,  w)  Plsin<A 
ij(&>)|,  |£(<u):  • » |^(„)  C(oj)lpf 


Lim  w(x. 
jq(a>)j,  K(o>)  • 


-BlPlx,  <,j)  P|cos<i 


(26) 

(27) 


Comparison  of  genet  al  calculations  with  observations 

To  test  eq  26  and  27.  we  need  to  determine  the  coherence  (and  phase  lag)  at  different  fre- 
quencies between  the  ice  deformation  time  series  and  the  atmospheric  pressure  time  series.  In 
particular,  we  would  like  to  estimate  ijL.d  and  £(i.»).  To  carry  out  such  an  estimation,  we  note  that 
for  a linear  system  the  frequency  response  function  may  be  estimated  by  a cross  spectral  analysis 
(Jenkins  and  Watts  1968,  p.  352).  Using  the  unfiltered  time  series  \<t),  <o(t)  and  Pit)  (the  camp 
atmospheric  pressure),  a cross  spectral  analysis  was  carried  out  using  the  lagged  product  method. 
Figure  4 shows  the  resulting  coherence  spectra  and  phase  angles.  The  phase  angle  convention  is 
such  that  a positive  phase  angle  indicates  a deformation  signal  lagging  behind  the  atmospheric 
pressure.  Using  eq  26  and  27  as  models,  a negative  phase  angle  of  90  would  occur  for  a perfect 
Hooke’s  law  behavior.  Figure  5 illustrates  the  resulting  amplitudes  of  //(<,»)  and  £(*>)  obtained  from 
estimates  of  the  amplitude  of  the  response  functions  of  A and  >.i  upon  P. 

Figure  5 shows  that  both  the  bulk  and  shear  "viscosities"  £(c»)  and  >/Ln)  exhibit  a general 
decrease  in  amplitude  with  increasing  frequency.  Referring  to  Figure  4.  the  phase  angle  behavior 
does  show  some  negative  tendency  indicative  of  elasticity,  especially  in  the  vorticity-pressure 
phase  at  higher  frequencies.  However,  the  overall  behavior  would  generally  seem  to  be  more 
suggestive  of  a viscous  behavior  (0  phase  angle)  rather  than  an  elastic  ( 90  phase  angle)  be- 
havior. 

The  decrease  in  the  "viscosity”  amplitudes  with  increasing  frequency  is  plausible  on  physical 
grounds;  for  example,  let  us  imagine  farcing  a simultaneous  sinusoidal  oscillation  in  the  divergence 
rate  and  shear  rate  of  a given  region  of  pack  ice.  The  displacements  of  the  oscillations  will  scale 
as  1/oj.  Consequently,  at  very  low  frequencies  the  average  compressive  stress  magnitude  over  one 
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Figure  4.  Coherence  spectra  and 
phase  angles  between  a)  vorlicity  and 
atmospheric  pressure  and  b)  divergence 
rate  and  atmospheric  pressure.  The 
95%  conlidence  limits  lor  the  phase 
angles  vary  from  t 20°  to  * 25°. 
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Figure  5.  Frequency  dependence  ol  the 
generalized  bulk  and  shear  viscosity  am- 
plitudes. The  amplitudes  were  obtained 
using  the  estimated  response  function  ol 
A and  w upon  P. 


cycle  should  be  larger  than  that  at  higher  frequencies  because  the  larger  compressive  displace- 
ments might  cause  more  thick  ice  to  be  crushed.  Similarly,  the  shear  stress  should  be  larger  for 
highly  compressed  ice;  and  since  it  is  observed  that  shear  and  dilatation  deformations  generally 
occur  simultaneously  (Ilibler  et  al.  1974),  the  average  shear  stress  magnitude  for  a cycle  might  be 
expected  to  increase  with  .he  increasing  displacement  amplitude  of  the  lower  frequency  cycles. 

Referring  back  to  Figure  3,  we  see  that  most  of  the  atmospheric  pressure  variance  is  at  wave- 
lengths longer  than  100  hours  (90%  of  the  pressure  variance  in  Fig.  3 is  at  wavelengths  greater 
than  120  hours).  The  fact  that  £(«),  and  especially  are  relatively  flat  over  these  wavelengths 
(factor  of  3 variation),  coupled  with  the  small  phase  shifts  at  low  frequency,  indicates  that  at  low 
frequencies  the  Glen  law  limiting  case  li.e.,  i;b,»)  ana  £(a»)  constant]  is  a good  first  approximation 
to  the  generalized  constitutive  law.  Since  a linear  law  of  any  kind  is  probably  only  a crude  approxi- 
mation of  the  true  ice  rheology,  the  above  results  suggest  that  for  predictive  purposes  there  is 
little  advantage  in  using  a more  generalized  linear  law  as  opposed  to  the  Glen  law. 


Conclusions 

The  most  obvious  inadequacies  of  the  comparison  made  here  are  the  neglect  of  finite  boundaries 
in  the  drift  predictions  and  the  use  of  a simplified  ice  rheology.  However,  the  calculations  and 
comparisons  in  this  paper  do  provide  some  helpful  insight  into  expected  differential  sea  ice  drift  for 
different  ice  conditions.  Specifically,  this  study  indicates  several  conclusions  relative  to  AIDJEX. 
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1.  The  general  agreement  between  the  infinite  boundary  linear  drift  theory  predictions  and 
observations  indicates  that  the  dominant  aspects  of  the  mesoscale  differential  drift  behavior 
observed  in  the  1971  and  1972  A1DJEX  pilot  programs  may  lie  explained  using  simple  boundary 
conditions  and  straightforward  linear  constitutive  laws.  Certainly  more  complete  calculations  are 
needed  to  explain  detailed  drift  behavior. 

2.  The  solution  of  the  linear  drift  calculation  indicates  the  sensitive  nature  of  differential 
comparisons  in  that  smaller  values  of  the  constitutive  law  parameters  r/  and  C do  not  only  change 
the  magnitude  of  the  divergence  rate  but  completely  change  its  sign. 

3.  With  respect  to  spatial  scales,  the  infinite  boundary  linear  drift  solution  indicates  that 
for  long  wavelength  variations  in  ice  deformation  the  internal  ice  stress  is  unimportant,  whereas 
for  short  wavelength  variations  the  internal  ice  stress  becomes  critical.  Consequently,  ic  drift 
calculations  using  the  same  viscous  parameters  may  be  quite  different  for  pressure  fields  varying 
slowly  in  spree  as  opposed  to  those  varying  rapidly.  For  example,  given  appropriate  ./  and  l, 
values,  it  is  possible  to  have  the  ice  converge  in  a high  pressure  system  covering  most  of  the 
Arctic  Basin  and  to  diverge  in  a high  pressure  system  covering  only  a portion  of  the  Arctic  Basin. 

Such  spatial  scaling  effects  may  also  have  a bearing  on  comparisons  of  different  time  scales, 
since  the  temporal  mean  of  the  pressure  field  over,  say  a month,  may  vary  more  slowly  in  space 
than  the  mean  daily  pressure. 

1.  The  cross  spectral  study  of  the  atmosphere  pressure  and  ice  deformation  in  light  of  a 
generalized  linear  constitutive  law  indicates  that  our  observed  strain  results  may  be  better  ex- 
plained by  a viscoelastic  law  including  memory  effects  than  by  a simple  viscous  linear  law.  How- 
ever, at  low  temporal  frequencies  (<  0.01  hr-1)  the  generalized  Hooke’s  law  is  similar  to  the  Glen 

viscous  law,  this  suggests  that  for  general  predictive  purposes  such  a generalized  linear  law  will 
probably  yield  only  slight  improvement  over  the  Glen  law. 

■r).  Finally,  the  fact  that  differential  drift  follows  the  local  pressure  field  reasonably  well 
indicates  that  the  ice  velocity  field  may  be  quite  nonlinear.  Consequently,  differential  drift  esti- 
mation using  long  strain  lines  ( 100  km)  may  not  adequately  resolve  high  wave  number  variations 

in  the  ice  velocity  field. 
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APPENDIX.  RELATIVE  MAGNITUDES  OF  DIFFERENTIAL  DRIFT  FORCES 

A substitution  of  measured  drift  parameters  into  eq  6 and  7 tf  Part  11  of  this  report)  allows  a 
direct  assessment  to  be  made  of  the  relative  magnitudes  of  the  wind,  water  and  Cmohs  Mrt 
terms  independent  of  the  value  and  functional  form  of  the  internal  ice  stress  and  nidt  pendent  o 

boundary  conditions.  Figure  2 shows  that  typical  values  for  A * and  V • v wind  ^e  given  by 
\ 0 0002  H1  1 w 0.0006  hr-4.  V . y 0.14  hr-1  . Using  eq  6 and  values  of  B,  D,  0,  <i>  and  f 

as  mentioned  earlier,  we  find  that  the  wind  stress  term  is  about  10  to  20  times  as  large  as  the 

water  and  Coriolis  stress  terms.  This  indicates  that  for  dlff“ent^ 

strain  and  Coriolis  terms  for  compact  conditions  is  reasonable,  and  thus  jnstif.es  the 

large  r/,£  drift  solution. 

1,  is  useful  to  contrast  the  differential  drift  results  in  this  paper  with  regular  drift  results 
For  regular  (nondifferent ial)  drift,  the  wind  to  ice  velocity  ratio  is  typically  of  the  order  of  o0  or 
less  (Reed  and  Campbell  1962,  Skiles  1968).  Thus,  since  D B 50.  for  regular  drift,  water  an 
Coriolis  stress  terms  may  not  be  neglected.  For  differential  drift  results  in  this ^paper  on  t 
other  hand,  the  ratio  of  wind  divergence  to  ice  divergence  or  vorticitv  is  300.  so  that  water 
Coriolis  stress  are  relatively  small. 
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